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In order to derive equation of state for the pure 5(7(3) Yang-Mills fields from first principles, 
it is proposed to generalize the effective potential approach for composite operators to non-zero 
temperatures. It is essentially non-perturbative by construction, since it assumes the summation 
of an infinite number of the corresponding contributions. There is no dependence on the coupling 
constant, only a dependence on the mass gap, which is responsible for the large-scale structure of the 
QCD ground state. The equation of state generalizes the Bag constant at non-zero temperatures, 
while its nontrivial Yang-Mills part has been approximated by the generalization of the free gluon 
propagator to non-zero temperatures, as a first necessary step. Even in this case we were able to 
show explicitly that the pressure may almost continuously change its regime at T* — 266.5 MeV. All 
the other thermodynamical quantities such as energy density, entropy, etc. are to be understood to 
have drastic changes in their regimes in the close vicinity of T* . All this is in qualitative and quan- 
titative agreement with thermal lattice QCD results for the pure Yang-Mills fields. We have firmly 
established the behavior of all the thermodynamical quantities in the region of low temperatures, 
where thermal lattice QCD calculations suffer from big uncertainties. 

PACS numbers: 11.15.Tk, 12.38.Lg 



I. INTRODUCTION 

The prediction of a possible existence of the Quark-Gluon Plasma (QGP) was one of the best theoretical achieve- 
ments of Quantum Chromodynamics (QCD) at non-zero temperatures and densities (the rather full list of the cor- 
responding references can be found in the text-book on finite-temperature field theory in Ref. [l| and in Ref. as 
well). The equation of state (EoS) for the QGP has been derived analytically up to the order g^\n{l/g'^) by using 
the perturbation theory (PT) expansion for the evaluation of the corresponding thermodynamical potential term by 
term ([l|,[3,[l] and references therein). 

However, the most characteristic feature of this PT expansion is its non-analytical dependence on the QCD coupling 
constant g^. In fact, this means that the PT QCD is not applicable at finite temperatures, apart from maybe at very 
high temperatures. The problem is not in poor convergence of this series [H, [IjSl (^"^ mathematics there exist methods 
how to improve convergence). The problem is that in the case of the above-mentioned non-analytical dependence 
one cannot even define the radius of convergence, so any next calculated term can be bigger than the previous one. 
This is the principle problem which can be resolved by no means. From the strictly mathematical point of view 
the four-dimensional QCD at non-zero temperatures effectively becomes the three-dimensional theory. At the same 
time, the three-dimensional QCD has much more severe infrared singularities [4] and its coupling constant becomes 
dimensional. That is the reason why the dependence becomes not analytical, while using the dimensionless coupling 
constat g^ (one needs to introduce three different scales, T, gT and g^T, where T is the temperature in order to 
somehow understand the dynamics of the QGP within the thermal PT QCD approach). Thus there is an exact 
indication that the analytical EoS derived by thermal PT QCD is wrong. 

At present, the only method to be used in order to investigate thermal QCD is the lattice QCD at finite temperature 
and baryon density which underwent a rapid recent progress (P, [3, S S 0| and references therein). However, the 
lattice QCD, being a very specific regularization scheme, first of all is aimed at obtaining the well-defined corresponding 
expressions in order to get correct numbers from them. So, one gets numbers, but not understanding on what is going 
on. Such kind of understanding can only come from the dynamical theory which is continuous QCD. For example, 
any description of the QGP is to be formulated in the framework of the dynamical theory. The lattice thermal QCD is 
useless in this. The need in the analytical EoS remains, but, of course it should be essentially non-perturbative (NP), 
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FIG. 1: The skeleton loop expansion for the efltective potential. The wavy lines describe the full gluon propagators D, while the 
solid lines - the full quark propagators S. T is the full quark-gluon vertex, while T3 and T4 are the full three- and four-gluon 
vertices, respectively. The ghost skeleton loops are not shown explicitly. 



reproducing the thermal PT QCD results at a very high temperature only. Thus analytic NP QCD and lattice QCD 
approaches to finite-temperature QCD do not exclude each other, but contrary they should complement each other. 



results of SU{i) lattice QCD calculations for the QGP EoS. 

The formalism we are going to use in order to generalize it to non-zero temperature is the effective potential approach 
for composite operators [9| . It is essentially NP from the very beginning, since it is dealing with the expansion of the 
corresponding skeleton loop contributions (for more detail description see section 2 and Ref. as well, where it has 
been generalized on quark degrees of freedom, but not using the confinement-type solution for the quark propagator). 
The main purpose of this paper is to derive EoS for the gluon matter by introducing the temperature dependence 
into the effective potential approach in a self-consistent way, in particular by using the confinement-type solution for 
the full gluon propagator (see below). 



The quantum part of the vacuum energy density (VED) is determined by the effective potential approach for 
composite operators [§]. In the absence of external sources the effective potential is nothing but the VED. It is given 
in the form of the skeleton loop expansion, containing all the types of the QCD full propagators and vertices, see Fig. 
1. So each vacuum skeleton loop itself is a sum of an infinite number of the corresponding PT vacuum loops, i.e., it 
contains the point-like vertices and free propagators (see Fig. 2, where one term only in each lower order is shown, 
for simplicity). The number of the vacuum skeleton loops is equal to the power of the Planck constant, fi. 

Here we are going to formulate a general method of numerical calculation of the quantum part of the truly NP 
Yang-Mills (YM) VED in the covariant gauge QCD. The gluon part of the VED to leading order (the so-called log- 
loop level ^ fi,, the first skeleton loop diagram in Fig. 1, and which PT expansion is shown explicitly in Fig. 2) is 
analytically given by the effective potential for composite operators as follows [9i]: 



where D{q) is the full gluon propagator and DQ{q) is its free counterpart (see below). Traces over space-time and 
color group indices are assumed. Evidently, the effective potential is normalized to V{Do) = 0. Next-to-leading and 
higher order contributions (two and more vacuum skeleton loops) are suppressed at least by one order of magnitude 
in powers of h. They reproduce very small numerical corrections to the log-loop terms, and thus are not important 
for the numerical calculation of the VED to leading order. 




II. THE VED 




(2.1) 




+ 



FIG. 2: Infinite series for the gluon part of the VED (the first skeleton diagram in Fig. 1) 
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The two-point Green's function, describing the full gluon propagator, is 



D^,{q) = ~i {T^Md{-q\0 + iL^M) ^, (2.2) 
where ^ is the gauge-fixing parameter and 

T^u{q) = gt,u -'^^ = g^^u- L^,^{q). (2.3) 

Its free PT counterpart Do = ^"ly('Z) is obtained by putting the full gluon form factor d{~q^ , in Eq. (2.2) simply 
to one, i.e., 

D%{q) = {T^u{q) + iL^u{q)} ^. (2.4) 
In order to evaluate the effective potential (2.1), on account of Eqs. (2.2) and (2.4), we use the well-known expression 

TrlniD^'^D) = 8 x AlndetiD^^ D) = 32 ln[(3/4)d(-g2, + (1/4)]. (2.5) 

It becomes zero indeed when equating (i(— g^,^) — 1. 

Going over to four-dimensional Euclidean space in Eq. (2.1), one obtains (e^ — V{D)) 



= -16 



(2^)4 



In[l + 3d(g2,0]-^d(g2,0+a 



(2.6) 



where the constant a = (3/4) — 21n2 — —0.6363 and the integration from zero to infinity is assumed. The VED Eg 
derived in Eq. (2.6) is already a colorless quantity, since it has been already summed over color indices. Also, only 
the transversal ("physical") degrees of freedom of gauge bosons contribute to this equation, so there is no need for 
ghosts to cancel their longitudinal (unphysical) counterparts. 

However, the derived expression (2.6) remains rather formal, since it suffers from different types of the PT contri- 
butions ("contaminations"). In order to define the truly NP VED free of all the above-mentioned problems, let us 
make first the identical transformation of the full effective charge in Eq. (2.6) as follows: 



d{q\0 = d{q\0 - dP^{q\0+dP^{q^^) = ^ {q^) + d^^ {q\ £) . (2.7) 

where d^^{q'^, ^) correctly describes the PT structure of the full effective charge d{q'^, ^), including its behavior in the 
ultra-violet (UV) limit, compatible with asymptotic freedom (AF) phenomenon in QCD [ij]), otherwise remaining 
arbitrary. On the other hand, d^^{q^) defined by the above-made subtraction, is assumed to reproduce correctly the 
NP structure of the full effective charge, including its asymptotic in the deep infrared (IR) limit. This underlines 
the strong intrinsic influence of the IR properties of the theory on its NP dynamics. Evidently, both terms are valid 
in the whole energy/momentum range, i.e., they are not asymptotics. Let us also emphasize the principle difference 
between d{q^,£) and d^^{q^). The former is NP quantity "contaminated" by the PT contributions, while the latter 
one, being also NP, is, nevertheless, free of them. Thus the separation between the truly NP effective charge d^^{q^) 
and its nontrivial PT counterpart d^^{q^,£,) is achieved. For example, if the full effective charge explicitly depends 
on the scale responsible for the truly NP dynamics in QCD, say A^p, then one can define the subtraction 



d^'-iq^Alp) = d{q\Klp) - d{q\Klp = 0) = d{q\Klp) - d^^(g2), (2.8) 

which is obviously equivalent to the decomposition (2.7). In this way the above-mentioned separation becomes exact 
and unique as well (for such concrete example see below). Let us emphasize that the dependence of the full effective 
charge d{q^ , ^%ip) ^a^p ^^^^ only regular. Otherwise it is impossible to assign to it the above-mentioned physical 
meaning, since A^p can be only zero (the formal PT limit) or finite, i.e., it cannot be infinitely large. In principle, 
in some special models of the QCD vacuum, for example such as the Abelian Higgs model [H, [ij] , the NP scale is to 
be identified with the mass of the dual gauge boson. Let us note that if there is no exact criterion how to distinguish 
between the truly NP and the nontrivial PT parts in the full effective charge as described above, then it is possible 
from the full effective charge to subtract its UV asymptotic only. Evidently, in this case the separation between the 
truly NP and the nontrivial PT parts may not be unique. 
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III. GENERALIZATION TO NON-ZERO TEMPERATURES 

Substituting the above-discussed exact decomposition (2.7) into Eq. (2.6), introducing further the effective scale 
squared, separating the NP region from the PT one (soft momenta from hard momenta), and omitting some algebraic 
rearrangements (see Refs. [l^. fl5| and especially recent paper for details), one obtains 



eyM(r) = -Bym + Bym{T) + Pym{T). 

Here evidently Cg = eyM and Bym is the Bag constant at zero temperature [l( 
explicitly given by the following expressions 



(3.1) 

Also, Bym{T) and Pym{T) are 



SyM(T) = 16 



=// (Tq 



ln[l + 3ar(a-^ar(<Z^) 



(3.2) 



and Pym{T) has more complicate form, namely 



P.M(T)^ -16/^ 



d^ 



ln[l + 3ar(g2) + ia^m " IW.^ i^) + aril')] + « 



respectively, since it depends on both effective charges. In all these equations 



jNP 



,PT 



iPT 



(1% 



(3.3) 



(3.4) 



because d^^{q'^) and d^^{q^) are the truly NP and the nontrivial PT effective charges, respectively, as it follows from 
above. Precisely these expressions should be generalized to non-zero temperatures in order to get EoS for the pure YM 
fields. That is why we introduce the dependence on the temperature T in advance. Evidently, Eq. (3.2) will reproduce 
the temperature-dependent Bag constant. In the expression for Pym{T) the integration is from zero to infinity, while 
in the integral for Bym{T) it is from zero to the effective scale squared q^f f, which just symbolically shown in Eq. (3.2). 
It is worth emphasizing that a so defined Bag constant (3.2) is free of all types of PT contributions (" contaminations" ) , 
as it is required (this was a reason for the above-mentioned algebraic rearrangements and subtractions, see Ref. [l6| 
and references therein). 

The problem remaining to solve is to choose the truly NP effective charge aril')- For the different truly NP 
effective charges we will get different analytical and numerical results. That is why the choice for its explicit expression 
should be physically and mathematically well justified. Let us choose the truly NP effective charge as follows: 



NP 



is? 



(3.5) 



where K^p is the mass scale parameter (the mass gap) responsible for the large-scale structure of the true QCD 
vacuum. It is well known that in continuous QCD it leads to the linear rising potential between heavy quarks, " seen" 
by lattice QCD [l3, [l3| as well (('7^)~^-type behavior for the full gluon propagator). Moreover, in Ref. [l0| it has 
been explicitly shown that it is a direct nonlinear iteration solution of the transcendental equation for the full gluon 
propagator in the presence of a renormalized mass gap. The separation between the truly NP and the nontrivial PT 
effective charges is both exact and unique, since the PT effective charge is always regular at zero, while the truly NP 
effective is singular at the origin (in the formal PT limit (A^p 0) the truly NP effective charge vanishes, while 
its nontrivial PT counterpart will survive). Let us also note that the chosen effective charge (3.5) does not depend 
explicitly on the gauge choice. It has been already used [ij, [3, [j^ in order to calculate the Bag constant, which 
turned out to be in a very good agreement with such important phenomenological parameter as the gluon condensate. 
It leads to many other desirable properties for the Bag pressure at zero temperature [lB|. Thus, our choice (3.5) is 
physically justified and mathematically confirmed, as required above. 

In the imaginary time formalism [J [2^ , these expressions can be easily generalized to non-zero temperatures T 
according to the prescription (let us remind that there is already Euclidean signature) 



dgo 
(2^) 



T 



+00 

n— — oo 



9^ = q^ 



ql = Ci^ + ujI = uj"^ + ul, 



ujn = 2m:T, 



(3.6) 
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i.e., each integral over go of the loop momentum is to be replaced by the sum over Matsubara frequencies labelled by 
n, which obviously assumes the replacement qo ^ ojn = 2mrT for bosons (gluons). In frequency-momentum space 
the truly NP effective charge becomes 



A2 

q2 +a;2 



PI r 2\ PI / 2 2\ PI I 2 2\ 



It is also convenient to introduce the following notations 



(3.7) 



(3.8) 



where, evidently, in all expressions here and below is the three-dimensional loop momentum squared in complete 
agreement with the relations (3.6). 

IV. THE DERIVATION OF Bym{T) 
In frequency-momentum space the Bag pressure (3.2) after the substitution of the relations (3.6) becomes 



n— — oo 



ln[l + 3af-(q^a;^)]-^af-(q^a;^) 



(4.1) 



where the truly NP effective charge is given in Eq. (3.7), and notations of Eqs. (3.7)-(3.8) are also valid, of course. 
After its substitution into Eq. (4.1), one yields 



/i3 j-oo 
* ' Ti — — rx- 



H3A%P + q2 + col] - ln[q' + ^l] - l^NP ^l^o 



(4.2) 



The summation over the Matsubara frequencies squared — [IttTYh? can be easily done, and the dependence on 
the effective scale oj^f f (see Appendix) is omitted, for simplicity. Here it is also convenient to introduce the following 
notation 



CO' = ^q2 + m'^^f = ^q2 + ^K^^ = ^u^^ + 3A^p = c.^1 + ^A^ 



'■NP-: 



(4.3) 



So it is possible to say that within our approach to non-zero temperatures we have two sorts of gluons: massless uo 
and massive uo' with the effective mass 



In the second term the summation over Matsubara frequencies can be done explicitly, namely 



(4.4) 



+ 00 -| oo 

V t = 



1 



+ 00 



q" + < 



n— — oo 



a;2 (27rT)2n2 



= (2,r//3)-2 



1 



n— — OD 



n2 (/3a;/27r)2 



= (27r//3)-'(27rV/3a;) 1 + 



Ari + ^- 
2a; I e^'^ -\ 



(4.5) 



A. The summation of logarithms 



In terms of the above-introduced parameters the sums in Eq. (4.2), containing the corresponding logarithn 
like 



+ 00 



and 



+00 



ln[q' + ujI] = lna;2 + 3 ^ ln{27,/pf[n^ + (/3a;/27r)2]. 



It is convenient to introduce the notations as follows: 



L{tj') = J2 Hn^ + Wl'^T^f] = ^ Inn^ + ^ In 



n=l n=l 



and equivalently 



L{uj) = Hn^ + (/?w/27r)2] = ^ Inn^ + ^ In 



71—1 n—1 



1 - 



Calculating explicitly the first one, we can calculate automatically the second by simply replacing oj' 
vice-versa. Evidently, in these expressions we introduce the following notations: 



PLU' 



~2 



So the difference L{u') — L{u) becomes 



L{uj') - L{uj) = ^ In 
or equivalently 



n=l 



El- 

n=l 



Insina;' — Ina;'^ — Insina; -|- ^ Ina;^, 



i(a;0-L(a;) = -^lnfe)+ln 



From the relation (4.10) it follows 



smx 
sin a; 



so the previous equation (4.12) finally becomes 



1 1 11- e~'^'^' 

Liu;') - L(a;) = - - In - + -/3(u;' - a;) + In ^--^ 
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B. The explicit expressions for the integrals 

Substituting all our results of the summations into Eq. (4.2), dropping a /3-independent terms plj], and performing 
almost trivial integration over angular variables, one obtains 



Bym{T) = — 



duj lJ^ 



3,9 1 1 1 (\-e-<^'^' 

Alp ^^-2/3-1 In' 



4 "'^ ^^^\l-e~P^ 

It is convenient to present the integral (4.15) as a sum of a few terms 



(4.15) 



ByM(T) = ~^^klpB%{T) \b%{T) ~ bP,,{T) 



where the explicit expressions of all these integrals are given below 



(4.16) 



duj- 



(4.17) 



(4.18) 



b{?]^(T) = '" dw In (l - e-^"') . 
In all these integrals the upper limit uj^ff is explicitly shown now and (3~^ = T, while 



(4.19) 



UJ ^ \ LU 



;2 + 3A^P = c.\/l + ^. 



(4.20) 



V. THE DERIVATION OF Pym{T) 

The term which contains the information about the nontrivial YM part (3.3) of the future gluon plasma EoS is 



Pym{T) = -16 



(2^)4 



ln[l + {l') + 3ar (9^)] - iK' (q') + «f + 



(5.1) 



where af^{q^) is the nontrivial PT effective charge. Due to the above-mentioned normalization of the effective 
potential approach in Eq. (2.1), the investigation of this part makes sense to begin with the approximation of the 
nontrivial PT part by its free PT counterpart, i.e., to put af'^{q^) = 1, as a first necessary step. Then Eq. (5.1) for 
the YM pressure in frequency-momentum space becomes 



+ 00 r 



/r] T 

and after substituting of the relations (3.7) into it, one obtains 



ln[l + Jar(q^c.,^J]-Jar(q^-,^J 



(5.2) 



Pym(T) = -16 



/,3 +00 
A^T y 



(5.3) 
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and the summation over the Matsubara frequencies squared = (27rT)^n^ can be easily done. For this purpose it 
is also convenient to introduce the following notation 



q2 + fhlf^ = \l<f + 1^%P = + 1^%P = 1 + ^a^^P' (5-4) 



so it is possible to say that within our approach to non-zero temperatures at this intermediate stage we have two sorts 
of gluons: massless w and massive cj with the effective mass 

rneff = ^J^NP = ^rn'^ff (5.5) 

Comparing Eqs.(4.2) and (5.3) one can write down the final result directly. For this purpose, in the final system 
of Eqs. (4.16)-(4.19) one must change the overall sign, replace u' by u> and integrate from zero to infinity. Thus, one 
obtains 



6 ,2 — 



Pym{T) = ^A^pP^^^(r) + [P^liT) - 41(T)J , (5.6) 
where the explicit expressions of all these integrals are given below 

POO 

P^L{T)=l dc.^^, (5.7) 

pPm{T) = dwu^ In (1 - e-'^-) , (5.8) 
Jo 

POO 

P^%{T)= / da;a;2ln(l-e-^^). (5.9) 
Jo 

VI. THE GLUON MATTER EOS 

Denoting further cyMiT) + Bym = Pgm{T) in the left-hand-side of our EoS (3.1), one obtains 

Pgm(T) = Bym{T) + Pym{T), (6.1) 

and in this equation Bym{T) and Pym{T) are given in Eqs. (4.16) and (5.6), respectively. Summing up all the 
integrals (4.17)-(4.19) and (5.7)-(5.9), one obtains that the gluon matter EoS (6.1) finally becomes, 

fi 1 fi 

Pgm{T) = ^A%pP,{T) + -iT[P2{T) + Ps{T) - P^T)], (6.2) 
where the dependence on the thermodynamical variable T is only shown explicitly and 

P,{T) = r du;^^ (6.3) 



while 



P3(r) = ^ °" doj uj^ In (l - e"'^'^') 

/■OO 

P^T) = / da;a;2ln(l-e-''^). 
Jo 



(6.4) 
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Let us recall once more that in all integrals f3 = and ujeff is fixed (see Appendix), while 



Oj^^uj^ + ^A%p, co' = ^uj^ + 3A%p. (6.5) 

In the formal PT limit (A^p — 0) from these relations it follows that Hi ^ oj' = uj and the combination P2{T) + 
PsiT) — Pi{T) becomes identical zero. Thus the gluon matter pressure (6.2) in this limit vanishes, i.e., it is truly NP, 
indeed. 

The effective potential has been normalized to zero in the D — )■ Dq limit, which reproduces the so-called Stefan- 
Boltzmann (SB) non-interacting (ideal) gas of massless particles (gluons) at high temperatures So the SB limit 
can be added (if necessary) to the truly NP pressure (6.2) in the T ^ oo limit only, i.e., 

Pgm{T) ^ PsB = ^T^^T\ T ^ oo (/? 0). (6.6) 
45 

In the same way, the corresponding SB limits should be added (if necessary) to all other thermodynamical quantities 
considered below. 



VII. THERMODYNAMICAL POTENTIAL AND OTHER THERMODYNAMICAL QUANTITIES 



In quantum statistics the thermodynamical potential ^1{T) is nothing but the pressure P(T) apart from the sign, 
i,e, in our case we can put 



n{T) = ~Pgm{T). (7.1) 

In quantum statistical theory all the important quantities such as energy density, entropy, etc., are to be expressed 
in terms of the thermodynamical potential. However, in the truly NP approach we cannot use the trivial relations 
between them which traced back to the PT even at non-zero temperatures. So the general formulae which to be used 
are [l| 



.(T) = ,7.2) 

for the pure YM fields, i.e., when the chemical potential is equal to zero. Evidently, here and everywhere below e and 
s are energy density and entropy, respectively, of the pure NP gluon matter. 

A. The energy density 

From Eqs. (7.1)-(7.2) it follows that 

e{T) = Ts(T) - P(T), (7.3) 
so substituting the corresponding explicit expressions (7.1) and (6.2) and doing some algebra, one obtains 



e(T) = Aa^p {t^^ Pi(T)) + ^T'^Min (7.4) 

where obviously d = 9, since the only dependence on T is present. Here and everywhere below we introduced the 
following notation: 



M{T) = P2{T) + PsiT) - Pi{T). 



(7.5) 
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Also, here and below all the integrals and their derivatives can be explicitly obtained from the expressions (6.3)-(6.5). 
Again the SB energy density 

94 

esBiT) = -TT^T^ (7.6) 
should be added to our expression (7.4) in the high temperature T ^ oo (/3 0) limit only. 

B. The entropy 

In the same way the entropy (7.2) becomes 



6 dPi(r) 16 , , 16 d , , , , 

s{T) = -A^p— + -M{T) + -T—M{T), 7.7 

TT^ dJ TT^ TT^ dJ 



and again the SB entropy 

ssb{T) = I^^T^ (7.8) 
should be added to our expression (7.7) in the high temperature T ^ oo (/3 — > 0) limit only. 

C. The heat capacity 

One of the interesting thermodynamical characteristics of the QGP is the heat capacity cy, which is defined as the 
derivative of the energy density. Then from the thermodynamical relations (7.1)-(7.2) it follows 



dT \ dT ^ 

Using the explicit expression for the energy density (7.4), one finally obtains 



As in previous cases, the SB heat capacity 

4^{T) = I^^T^ (7.11) 
should be added to our expression (7.10) in the high temperature T — s- oo (/3 — > 0) limit only. 

VIII. NUMERICAL RESULTS AND DISCUSSION 

All our numerical results are present in Fig. 3. It is seen explicitly that the NP gluon pressure may almost 
continuously change its regime in the close neighborhood of a maximum at T* — 266.5 MeV in order to achieve 
the thermodynamical SB limit at high temperatures. For the displayed quantities in Fig. 3 the SB limits are the 
corresponding constants. At the same time, for all other thermodynamical quantities such as the energy density, 
entropy and heat capacity this is impossible as it follows from the curves shown in Fig. 3 (none of their power-type 
fall off at this point can be smoothly transformed into the constant behavior at high temperatures). In order to 
achieve the thermodynamical SB limits at high temperatures their full counter-parts should undergo drastic changes 
in their regimes in the close neighborhood of this point. As we already know from thermodynamics of SU{3) lattice 
QCD [Ei, [2H the energy and entropy densities have a discontinuity at a point Tc = 260 MeV, while the pressure 
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T, MeV 

FIG. 3: The NP pressure P, energy density e, entropy s and heat capacity cv as a functions of the temperature. The NP gluon 
pressure P has a maximum at T* = 266.5 MeV. 



remains continuous. Our characteristic temperature T* = 266.5 MeV is, surprisingly, very close to the same value. 
A clear evidence that something nontrivial in the behavior of the thermodynamical quantities in the vicinity of our 
characteristic temperature T* — 266.5 MeV should actually take place follows from the fact that at this point e — 3P, 
which should be valid at a very high temperatures only (SB limit). In other words, in order to derive EoS valid above 
T* , and thus to provide a correct picture of thermodynamics of the gluon matter in the whole range of temperature, 
one needs the nontrivial approximation of the YM part (5.1), compatible with AF phenomenon in QCD [11]. Just 
this will be subject of the subsequent paper. 

If we were not aware of the thermal lattice QCD results then we would be able to predict them. But we aware of 
them, so lattice results confirm our expectations of a sharp changes in the behavior of the entropy and energy densities 
in the region where the pressure is continuous. At the same time, it is worth emphasizing that we have no any problems 
in describing the behavior of all the important thermodynamical quantities at low temperatures below T* (see Fig. 
3). Moreover, apparently for the first time it is possible to predict their behavior in the region of low temperatures 
within our approach (there are no convincing lattice data for this region). We do not expect any serious changes in 
the behavior of the thermodynamical quantities in this region (exponential fall off or rise when the temperature goes 
down or up, respectively) even after taking into account the above-mentioned nontrivial approximation of the YM 
part (5.1), apart from "non-physical" maximums which should disappear, of course. However, whatever changes may 
occur they will be under our control. 

The confinement dynamics (3.5) generalized on zero-temperatures in Eq. (3.7) is still important especially in the 
region of low temperatures even up to the temperature at which all the important thermodynamical quantities may 
undergo drastic changes in their behavior (apart from the pressure). From the structure of our EoS (see Eqs. (6.2)- 
(6.5)), it clearly follows that below We// (which fixes T*, there is no explicit dependence between them, but rather a 
correspondence) and thus below T* we have mainly the massive gluon excitations uj' of the dynamical origin, which 
can be interpreted as the glueballs with masses m'^fj = V^Anp = 1.17 GeV. Above T* the gluon matter consists 
mainly of the free gluons to {uj gluons are artifacts due to the approximation of the nontrivial PT effective charge by 
its free PT counterpart in Eq. (5.1), as well as the above-mentioned " non- physical" maximums and hence their "tails" 
at high temperatures in Fig. 3). Just the confinement dynamics determines the phase transition from glueballs to 
"free" gluons and vice- versa in the case of SU{3) YM fields within our approach, indeed. 
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APPENDIX A: THE SCALE-SETTING SCHEME 

From the relations (3.6) it follows that in frequency-momentum space a possible free parameter of our approach is 
the effective scale 

^eff = ^qlff~Ujl, (Al) 

where we introduced the constant Matsubara frequency lDc, which is always positive. So We// is always less or equal 
to qeff of the four-dimensional QCD, i.e., 

^eff < qeff- (A2) 

One then can conclude that qef / is a very good upper limit for uj^f / . In this connection, let us recall now that the Bag 
constant Bym at zero temperatures has been successfully calculated at a scale q^^^ = 1 GeV^, in fair agreement with 
other phenomenological quantities such as gluon condensate [l^ . So let us fixed the effective scale ojef / as follows: 

^eff = qeff = 1 GeV. (A3) 
The mass gap squared A^p calculated just at this scale is equal to [31 

A?^P = 0.4564 GcV^ (A4) 

Thus, we have no free parameters in our approach. The confinement dynamics is nontrivially taken into account 
directly through the mass gap, and not through the Bag constant itself. 
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